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ANALYSIS  OF  THE  NORTH  TRURO  SEA  CLUTTER  DATA 


ABSTRACT 


\ 


In  this  report  we  present  the  results  of  the  analysis  of  sea— clutter  data  collected  at 
North  Truro,  Massachusetts  by  the  MIT  Lincoln  Laboratory.  The  data  were  obtained  in 
conditions  typical  of  sea  state  2  at  X-,  S-,  L-bands,  UHF  and  VHF.  We  examined  the 
data  both  spectrally  and  statistically.  The  results  showed  that  sea  clutter  exhibited  very 
different  spectral  characteristics  at  higher  frequencies  compared  to  those  at  lower 
frequencies.  Experimental  sea-clutter  coefficients  as  a  function  of  grazing  angle  at  various 
frequency  bands  were  obtained  for  up— swell  and  cross-swell  conditions.  Most  sea-clutter 
models  use  the  radar  look  direction  relative  to  the  wind  direction  as  an  input  parameter 
(e.g.,  up-wind,  cross-wind,  etc).  Hence  a  direct  comparison  of  sea-clutter  magnitude  with 
those  calculated  from  models  under  different  wind  conditions  was  not  made.  Instead,  we 
showed  these  results  together  with  those  calculated  from  two  sea-clutter  models,  assuming 
up— wind  and  cross-wind  conditions.  The  data  were  quite  dissimilar  in  the 
low-grazing-angle  regions  to  those  predicted  by  the  two  models.  .Fi^th^  rpsparrh  on 
dutter  in  the  low-grazing— angle  regions  is  needed  to  verify  and  explain  the  observed 
differences.  We  found  the  sea-dutter  amplitude  statistics  fit  the  K-distribution  best  in 
the  important  low-probability-of-false-alarm  region.  The  distribution  approached  that  of 
a  Rayleigh  model  for  vertical  polarization  and  low-resolution  waveforms.  However,  there 
was  a  significant  proportion  of  resolution  cells  with  sea  dutter  having  non-Rayleigh 
statistics  even  in  relatively  low-resolution  waveforms  and  vertical  polarization.  For 
horizontal  polarization  and/or  high— resolution  waveforms,  the  statistics  approached  those 
of  a  lognormal  model.  The  statistical  results  provided  further  evidence  that  the 
K-distribution  can  serve  as  a  limiting  distribution  for  sea  clutter.  Spatial  and  temporal 
characteristics  of  sea  spikes  were  examined  through  the  use  of  range-time-intensity  plots 
of  high-resolution  sea  dutter.  The  results  supported  the  theory  that  sea  spikes  are 
associated  with  the  wave  crests  in  sea  swells.  This  association  was  more  evident  from 
high-resolution,  V— POL  sea-clutter  data.  The  results  of  this  report  are  pertinent  to 
shipboard  radar  applications  where  the  grazing  angle  is  low.  They  also  provide  additional 
information  for  sea-clutter  modeling  in  the  low— angle  region,  particularly  at  the  lower 
frequency  bands. 
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ANALYSE  DES  DONNIES  DU  FOUILLIS  DE  MER  DE  NORTH  TRURO 


r£sum£ 


Dans  ce  rapport,  nous  presentons  les  resultats  de  l’analyse  des  donnees  du  fouillis  de 
mer  recueillies  a  North  Truro  au  Massachusetts  par  le  Laboratoire  Lincoln  du  MIT.  Les 
donnees  ont  £te  obtenues  dans  des  conditions  typiques  k  Petat  de  mer  2  et  aux  bandes  de 
frequences  X,  S,  L,  UHF  et  VHF.  Nous  en  avons  fait  Panalyse  spectrale  et  statistique.  Les 
resultats  ont  montre  que  les  caract4ristiques  spectrales  du  fouillis  de  mer  aux  hautes 
frequences  etaient  tres  differentes  de  celles  aux  frequences  plus  basses.  Des  coefficients 
experimentaux  pour  le  fouillis  de  mer  en  fonction  de  l’angle  d’incidence  rasante  a 
differentes  bandes  de  frequences  ont  ete  obtenus  pour  une  observation  face  aux  vagues  et  en 
travers  des  vagues.  Puisque  la  plupart  des  modeles  du  fouillis  de  mer  utilisent  comme 
parametre  d’entree  la  direction  d’observation  du  radar  par  rapport  a  celle  du  vent  (face  au 
vent,  en  travers  du  vent,  etc),  une  comparaison  directe  de  l’amplitude  du  fouillis  de  mer 
avec  celle  calculee  a  partir  de  ces  modeles  n’etait  pas  possible.  Nous  avons  plutot  compare 
ces  resultats  avec  ceux  calcules  a  partir  de  deux  modeles  de  fouillis  de  mer:  Pun  face  au 
vent  et  l’autre  en  travers  du  vent.  Les  donnees  correspondant  a  de  faibles  angles 
d’incidence  rasante  etaient  tr£s  differentes  de  celles  predites  par  les  deux  modeles.  Des 
recherches  supplement  air  es  sur  le  fouillis  de  mer  aux  incidences  rasantes  sont  necessaires 
pour  verifier  et  expliquer  les  differences  observees.  Nous  avons  trouve  que  la  distribution 
statistique  de  l’amplitude  du  fouillis  de  mer  ressemblait  le  plus  a  une  distribution  K  dans 
les  zones  a  faible  probability  de  fausses  alarmes  soit  celles  de  plus  grande  importance.  Pour 
des  formes  d’onde  a  faible  resolution  et  a  polarisation  verticale,  la  distribution  de 
l’amplitude  du  fouillis  s’approchait  de  la  distribution  de  Rayleigh.  Une  distribution  autre 
que  Rayleigh  a  ete  aussi  observee  pour  des  formes  d’onde  a  faible  resolution  et  a 
polarisation  verticale.  Pour  des  formes  d’onde  a  haute  resolution  ou  a  polarisation 
horizontale,  la  distribution  du  fouillis  s’approchait  d’une  distribution  log— normale.  Les 
donnees  statistiques  ont  fourni  des  preuves  supplementaires  que  la  distribution  K  peut 
servir  de  cas  Umite  pour  la  distribution  de  l’amplitude  du  fouillis  de  mer.  Les 
caracteristiques  spatio-temporelles  des  pointes  des  echos  de  mer  ont  ete  examinees  a  l’aide 
des  diagrammes  tridimentionnels  de  l’amplitude  du  fouillis  de  mer  en  fonction  de  la  portee 
et  du  temps  pour  des  formes  d’onde  4  haute  resolution.  Les  resultats  soutenaient  la  theorie 
selon  laquelle  ces  pointes  sont  associees  aux  cretes  des  vagues.  Cette  association  etait 
encore  plus  ^vidente  pour  des  formes  d’onde  a  haute  resolution  et  a  polarisation  verticale. 
Les  resultats  de  cet  article  sont  appropri£s  aux  applications  du  radar  maritime  ou  l’angle 
d’incidence  rasante  est  faible.  Ds  fournissent  egalement  de  Pinformation  supplemental 
sur  la  mod&isation  du  fouillis  de  mer  dans  les  zones  ou  l’angle  d’incidence  rasante  est 
faible,  surtout  aux  bandes  de  frequences  plus  basses. 
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EXECUTIVE  SUMMARY 


The  performance  of  shipborne  radars  is  adversely  affected  by  returns  from  the  sea 
surface  surrounding  the  ship.  These  returns  are  called  sea  clutter  in  radar  terminology. 
Although  sophisticated  signal  processing  techniques  have  been  developed  to  minimize  the 
interference  from  sea  clutter,  the  effectiveness  of  these  techniques  is  highly  dependent  on 
the  characteristics  of  sea  clutter.  In  order  to  optimize  shipboard  radar-signal— processing 
performance,  a  thorough  understanding  of  the  sea-clutter  behaviour  under  a  wide  variety 
of  sea  conditions  is  required.  There  are  two  important  characteristics  of  sea  clutter  which 
affect  radar  performance:  (a)  spectrum  and  (b)  temporal  amplitude-statistics.  These 
characteristics  change  as  the  radar  parameters  and  the  sea  conditions  change. 

We  present  in  this  report  the  results  of  the  analysis  of  sea-clutter  data  collected  at 
North  Truro,  Massachusetts  by  the  MIT  Lincoln  Laboratory.  The  data  were  obtained  in 
conditions  typical  of  sea  state  2  at  X-,  S-,  L-bands,  UHF  and  VHF.  We  examined  the 
data  both  spectrally  and  statistically. 

The  results  showed  that  sea  clutter  exhibited  very  different  spectral  characteristics 
at  higher  frequencies  compared  to  those  at  low  frequencies.  At  higher  frequency  bands 
(X— ,  S— bands),  the  sea  clutter  spectrum  was  generally  asymmetrical  with  a  longer  spectral 
tail  in  the  direction  of  the  Doppler  shift.  At  lower  frequency  bands  (L— band,  UHF  and 
VHF),  the  sea  clutter  spectrum  was  characterized  by  the  presence  of  double  spectral  peaks. 

Experimental  sea— clutter  coefficients  as  a  function  of  grazing  angle  at  various 
frequency  bands  were  obtained  for  up— swell  and  cross-swell  conditions.  Most  sea— clutter 
models  use  the  radar  look  direction  relative  to  the  wind  direction  as  an  input  parameter 
(e.g.,  up-wind,  cross-wind,  etc).  Hence  a  direct  comparison  of  sea-clutter  magnitude  with 
those  calculated  from  models  under  different  wind  conditions  was  not  made.  Instead,  we 
contrasted  these  results  with  those  calculated  from  two  sea-clutter  models,  assuming  an 
up-wind  and  a  cross-wind  conditions.  The  data  were  quite  dissimilar  in  the 
low-grazing— angle  regions  to  those  predicted  by  the  two  models.  Further  research  on  sea 
clutter  in  the  low-grazing— angle  regions  is  needed  to  verify  and  explain  the  observed 
differences. 

We  found  that  the  sea-clutter  amplitude  statistics  fit  the  K-distribution  best  in  the 
important  low— probability-of— false-alarm  region.  The  distribution  approached  that  of  a 
Rayleigh  model  for  vertical  polarization  ana  low-resolution  waveforms.  We  observed 
non— Rayleigh  statistics  even  in  relatively  low-resolution  and  vertical  polarization 
waveforms.  For  horizontal  polarization  and/or  high-resolution  waveforms,  the  statistics 
approached  those  of  a  lognormal  model.  The  statistical  results  provided  further  evidence 
that  the  K-distribution  can  serve  as  a  limiting  distribution  for  sea  clutter. 

Spatial  and  temporal  characteristics  of  sea  spikes  were  examined  through  the  use  of 
range-time-intensity  plots  of  high-resolution  sea  clutter.  The  results  supported  the 
theory  that  sea  spikes  are  associated  with  the  wave  crests  in  sea  swells.  This  association 
was  more  evident  from  high-resolution,  V— POL  sea-clutter  data.  This  result  also  offered 
some  clues  as  to  the  reason  why  the  Doppler  shifts  of  the  sea  clutter  appeared  to  be  greater 
for  horizontal  polarization  than  for  vertical  polarization,  at  least  at  X-  and  S-bands. 

The  results  of  this  report  are  pertinent  to  shipboard  radar  applications  where  the 
grazing  angle  is  low.  They  also  provide  additional  information  for  sea-clutter  modeling  in 
the  low-angle  region,  particularly  at  the  lower  frequency  bands. 
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1. 


INTRODUCTION. 


Sea  clutter  limits  the  performance  of  maritime  surveillance  radars.  Although  there 
exists  a  large  body  of  theoretical  and  experimental  work  in  the  literature  [1H4],  our 
understanding  of  sea-clutter  behaviour  is  still  incomplete.  A  thorough  understanding  of 
sea-clutter  behaviour  is  essential  if  optimum  performance  is  to  be  obtained  from  radar 
signal  processors.  For  example,  the  spectral  characteristics  will  influence  the  design  of 
moving  target  indicator  (MTI)  filters,  and  the  amplitude  statistics  of  sea  clutter  will 
influence  the  design  of  constant  false  alarm  rate  (CFARl  processors.  A  knowledge  of  the 
sea-clutter  strength  as  functions  of  various  system  ana  environmental  parameters  (e.g., 
frequency,  polarization,  waveform  resolution,  sea  state,  grazing  angle,  etc)  will  enable  the 
radar  designer  to  tailor  the  signal  processing  to  specific  conditions,  thereby  obtaining 
optimal  performance. 

In  this  report  we  present  the  results  of  an  analysis  of  the  sea-clutter  data  collected 
by  the  MIT  Lincoln  Laboratory  at  North  Truro,  Massachusetts.  This  set  of  data  was 
provided  to  Canada  under  a  Canada— U.S.  memorandum  of  understanding  in  clutter 
research.  The  results  provided  insights  into  the  various  mechanisms  which  contribute  to 
sea-clutter  returns  and  are  particularly  pertinent  to  shipboard  radars  where  the  grazing 
angle  is  low.  Section  2  describes  the  radar  equipment,  the  radar  site,  the  resulting  data 
base  including  data  calibration  and  the  signal  processing  techniques  employed  in  the 
analysis.  Section  3  describes  the  analysis  of  the  sea-clutter  data  and  discusses  the  results. 
Section  4  summarizes  the  findings  and  presents  concluding  remarks. 

2.  DATA  COLLECTION  AND  ANALYSIS. 

2.1  The  Radar  Equipment. 

The  sea— clutter  data  collection  was  carried  out  during  a  two— week  period  in  July 
1982  using  the  MIT  Lincoln  Laboratory  mobile  coherent  radar  facility [5].  The  technical 
characteristics  of  this  radar  are  summarized  in  Table  I. 


Table  I  :  Technical  Characteristics  of  the 
MIT  Lincoln  Laboratory  Mobile  Radar  Facility 

Frequency  (MHz)  165  430 

Azimuthal 

Beamwidth  (Deg)  13  5 

v - „ - / 

Range 

Resolution  36  m/150  m 

Polarization  Vert.  Transmit/Receive 
PRF  500  Hz 

RCS  Calibration  Accuracy  :  2  dB 


1250  3400  9100 

3  1  1 

N - v - - - ' 

15m  /  150m 
Hor.  Transmit/Receive 
to  4000  Hz 
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2.2  The  Radar  Site 


The  radar  site  was  located  at  the  North  Truro  Air  Force  Station  on  Cape  Cod, 
Massachusetts.  The  effective  radar  height  was  77.4  m  above  mean  sea  level.  An  ENDECO 
949  wave-track  buoy  was  leased  and  anchored  4  km  off  shore  at  40°  east  of  north  to 
provide  a  record  of  the  sea  condition  during  the  experiment  period.  The  wave  data  were 
sent  from  a  transmitter  on  board  the  buoy  to  a  shore-based  receiver  for  twenty  minutes  in 
each  hour.  These  data  were  processed  and  summarized  in  a  number  of  parameters  such  as 
significant  wave  height  (SWH)  and  significant  period  (SP). 

Figure  1  shows  the  significant  wave  height  during  the  experiment  period,  with  day 
”0"  being  July  2, 1982.  The  segments  in  the  figure  with  no  data  were  weekend  periods  when 
the  wave  buoy  receiver  was  not  in  operation.  The  ENDECO  949  did  not  measure  wave 
directions. 


DAY  OF  TEST 


Figure  1.  Variation  of  significant  wave  height  during  the  sea-clutter  experiment 
period. 


The  wind  velocity  readings  were  recorded  from  a  location  on  shore  at  the  site.  The 
averaged  wind  speed  was  about  5  knots  over  the  experimental  period  in  which  sea-clutter 
data  were  taken.  There  were  isolated  time  intervals  in  which  wind  speeds  in  the  8  to  13 
knots  range  were  registered. 


2 


2.3  The  Sea-Clutter  Data  Base 


The  significant  wave  height  for  most  experiments  ranged  from  0.3  m  to  0.45  m  (1  to 
1.5  ft).  Given  the  limitations  of  the  wave  data  and  wind  measurements,  we  could 
characterize  the  average  sea  condition  during  the  experimental  period  only  approximately, 
as  sea  state  (SS)  2  (with  reference  to  the  World  Meteorological  Organization  (WMO)  code 
[6]).  With  a  radar  neight  of  77.4  m,  the  grazing  angle  for  usable  data  was  below  5  degrees. 
Tne  experimental  facility  provided  for  (a)  coherent  data  collection,  (b)  multi-frequency 
measurement,  (c)  near-simultaneous  observation  of  the  sea  with  V—  and  H— polarizations 
and  (d)  long-dwell— mode  data  collection  permitting  spectral  and  statistical  measurements. 

There  were  approximately  350  data— collection  experiments,  defined  by  the  selection 
of  frequency,  polarization,  waveform  resolution  and  look  direction.  Data  were  taken  in  an 
azimuthal  sector  extended  from  0°  to  120°.  All  directions  were  measured  clockwise  from 
the  north  direction  (0°).  The  analysis  reported  here  was  carried  out  on  data  of  long-dwell 
experiments  that  employed  a  PRF  of  500  Hz  and  dwell  times  of  10  to  20  seconds.  Details 
of  this  data  base  may  be  found  in  [7]. 

Calibration 


MIT  Lincoln  Laboratory  personnel  recorded  the  coherent  sea-clutter  data  digitally 
and  calibrated  them  at  the  Laboratory.  Calibration  is  the  process  of  determining  the 
constant  which  relates  the  instantaneous  analog-to-digital  converter  (ADC)  readings  to 
the  effective  radar  cross-section  (RCS)  of  the  sea  clutter.  The  calibration  constant  was 
derived  from  the  two-way  radar  equation  using  in-situ  measurement  of  the  returns  from 
corner  reflectors  of  known  free-space  RCS  values.  The  corner  reflectors  were  mounted  on  a 
remote  tower  a  few  kilometers  from  the  radar.  These  calibration  experiments  were 
performed  for  each  frequency  band  and  polarization.  Since  propagation  of  radar  signals  on 
land  could  be  subjected  to  multipath  effects,  the  calibration  experiment  was  repeated  by 
translating  the  corner  reflector  vertically  along  the  length  of  the  remote  tower,  thereby 
tracing  out  the  interference  pattern.  The  appropriate  calibration  constant  was  then 
obtained  by  selecting  the  value  which  yielded  the  best  match  between  the  theoretical  and 
experimental  interference  patterns. 

2.4  Signal  Processing 

(a)  Spectral  Analysis 

We  computed  the  sea-clutter  spectra  using  the  modified  periodogram  method  [8]  for 
X — ,  S—  and  L-bands.  The  long-dwell  experiments  for  these  bands  employed  a  pulse 
repetition  frequency  (PRF)  of  500  Hz  and  a  dwell  length  of  10000  pulses.  This  provided 
nine  1024-point  subsequences.  We  applied  a  Hamming  window  to  each  subsequence  and 
performed  the  fast  Fourier  transform (FFT).  The  periodograms,  which  are  the  squared 
magnitudes  of  the  FFTs,  were  averaged  to  form  the  sea-clutter  spectrum. 

For  UHF  and  VHF,  we  used  a  specific  version  of  the  Maximum  Entropy  Spectral 
Estimation(MESE)  due  to  Burg  [9]to  obtain  an  improved  spectral  resolution.  The  length 
of  the  time  series  for  UHF  ana  VHF  was  5120  pulses  with  a  500  Hz  PRF.  At  UHF  and 
VHF,  the  spectral  width  of  sea  clutter  is  only  a  small  fraction  of  the  PRF.  The  spectral 
details  of  the  sea  clutter  at  these  frequency  could  be  lost  due  to  the  limited  resolution  of 
the  FFT  and  data  windowing  effect. 
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MESE  is  an  algorithm  which  provides  for  superior  spectral  resolution  for  random 
time  series  which  can  be  modeled  by  an  auto  regressive  (AR)  system.  The  details  of  MESE 
are  well  documented  £9]  JlO]  and  will  not  be  repeated  here.  For  the  results  presented  in 
Section  3,  the  MESE  (VHF)  spectra  were  obtained  using  an  AR  model  of  order  50. 

(b)  Decorrelation  Time 

The  decorrelation  time  of  the  sea  clutter  was  obtained  by  measuring  the  time  it  took 
the  envelope  of  the  sea-clutter  autocorrelation  function  to  decay  to  1/e  of  its  peak  value, 
where  e  is  the  base  of  natural  logarithm.  The  autocorrelation  function  of  a  complex 
sequence  {sn}  was  computed  by: 


ak  = 


N-l  s  s* 

£  sn  sn+k 
n=0 _ 

N-l  * 

E  s_  s_ 


(i) 


where  *  denotes  the  complex  conjugate. 


Since  there  could  be  ground-clutter  breakthrough  from  the  antenna  sidelobes,  the 
sample  mean  of  the  time  series  was  first  removed  before  correlation  was  performed.  This 
process  effectively  removed  any  ground-clutter  component  from  the  data. 

(c)  Estimation  of  Sea-Clutter  Coefficient 

The  sea-clutter  coefficient  of  a  resolution  cell  is  defined  as: 


<7°  =  a  /A 


(2) 


where  a  is  the  equivalent  radar  cross  section  of  the  sea-clutter  return  from  that  cell, 
and  A  is  the  area  of  the  resolution  cell  given  by 


(3) 


where  R  is  the  range,  $  is  the  azimuthal  beamwidth,  c  is  the  veto'  ty  of  light, 
r  is  the  pulse  length  and  <t>  is  the  depression  angle. 
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High  Cliitter-to-Noise  RatiofCNRl  Case 

For  data  with  high  CNR,  the  equivalent  sea-clutter  RCS  is  computed  by 
noncoherent  integration: 


1  N 

a  =  i  E  (Xi2  + 

N  i=l 


yi2) 


(4) 


where  Xi  and  yi  are  the  in— phase  and  quadrature-phase  samples  of  the  sea  clutter,  and  N  is 
the  number  of  samples. 

Low  CNR  Case 


From  spectral  analyses  of  the  North  Truro  sea-clutter  data,  we  found  that 
horizontally  polarized  sea  clutter  at  L— band  and  lower  was  orders  of  magnitude  lower  than 
that  for  vertical  polarization.  This  made  the  use  of  noncoherent  integration  for 
determining  clutter  coefficient  unacceptable  since  the  data  are  dominated  by  noise.  To 
remedy  this  problem,  a  for  the  low  CNR  case  was  determined  using  the  following  steps: 

(i)  Compute  the  sea-clutter  spectrum  from  the  data; 

(ii)  Determine  the  Doppler  interval  in  which  the  sea  clutter  spectrum 
lies(regions  where  sea-clutter  spectral  density  is  greater  than  the  noise 
density); 

(iii)  Obtain  a  by  summing  the  components  of  the  spectrum  over  the 

Doppler  interval  determined  in  (ii). 

The  elevation  beamwidths  of  the  X—  and  S— band  antenna  were  fairly  narrow  (about 
3.5°  for  X— band  and  4.7°  for  S— band,  one-way),  resulting  in  an  attenuation  of  returns  from 
close  ranges.  For  example,  at  a  range  of  2  km,  the  depression  angle  is  about  2.2°.  The 
S-band  elevation  gain  is  down  approximately  3  dB  at  this  angle.  As  a  result,  S-band  sea 
clutter  returns  at  this  range  will  suffer  an  attenuation  of  close  to  6  dB  (two-way).  In 
calculating  a  for  X—,  S— bands,  the  values  were  compensated  for  by  adding  the  amount  in 
dB  equal  to  the  negative  of  the  two-way  elevation  gain  of  the  antenna  at  the  depression 
angle  for  that  cell. 

3.  ANALYSIS  RESULTS 

3.1  Spectral  Analysis 

Figures  2a  and  2b  show  the  sea-clutter  spectra  for  two  X-band  experiments,  for 
vertical  and  horizontal  polarizations,  respectively.  Experimental  parameters  are  shown  on 
the  figures.  Values  of  SWH  and  SP  during  the  experiment  are  also  indicated.  The  spectral 
resolution  is  A  =  500/1024  =  0.488  Hertz  per  sample.  In  both  cases,  the  spectra  were 
obtained  at  a  range  ol  4150  m  in  three  look  directions  (30°,  75°  and  120°  azimuth).  The 
Doppler  shifts  of  the  sea-clutter  spectra  were  estimated  to  be  -40  Hz  (30°),  30  Hz(75°)  and 
50  Hz(120°)  for  the  V-POL  experiment  and  0  Hz(30°),  70  Hz(75°)  and  90  Hz(120°)  for  the 
H— POL  experiment. 
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SPECTRAL  DENSITY  (dBm*/AHz) 


DOPPLER  FREQUENCY  (Hz) 


(a)  V— POL 


DOPPLER  FREQUENCY  (Hz) 


(b)  H-POL 

Figure  2.  Typical  low -resolution  X-band  sea-clutter  spectra  in  the  up-swell  and 
near— cross— swell  directions. 
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The  spectral  power  was  substantially  greater  for  the  spectrum  of  the  120°  azimuth 
than  those  of  the  other  two  azimuths.  The  recorded  wind  speed  for  these  experiments  was 
7  knots  from  the  direction  of  220°.  This  did  not  correspond  to  the  direction  of  the 
maximum  Doppler  shift  in  the  sea-clutter  spectra.  There  was  no  consistent  correlation 
between  the  recorded  wind  direction  and  the  direction  of  maximum  Doppler  shift  in  the  sea 
clutter  throughout  the  experiment  period.  Since  ocean  current  is  known  to  be  present  in 
the  Cape  Cod  area,  it  may  be  inferred  that  the  observed  Doppler  shift  in  the  sea  dutter 
was  at  least  partly  the  result  of  ocean  current.  For  lack  of  a  more  appropriate  description, 
we  called  the  direction  with  the  maximum  positive  Doppler  shift  the  up-swell  (look 
direction  against  the  sea  wave)  direction  and  the  direction  with  near  zero  Doppler  shift  in 
the  sea-clutter  spectrum  the  cross-swdl  (look  direction  perpendicular  to  the  wave) 
direction.  This  result  indicated  that  the  up-swell /cross-swell  ratio  of  sea  clutter  was  high 
at  X-band. 

The  horizontal  polarization  experiment  (Figure  2b)  was  performed  about  5  minutes 
after  the  vertical  polarization  experiment.  Consequently,  the  sea  conditions  should  have 
been  quite  similar.  We  see  that  the  Doppler  shift  for  a  horizontal— polarization  sea-clutter 
spectrum  was  greater  than  that  for  vertical  polarization  for  similar  sea  conditions. 

Pidgeon  [11]  reports  this  phenomenon  and  gives  a  probable  explanation.  He 
attributes  the  difference  in  the  mean  Doppler  shift  between  the  two  polarizations  to  the 
different  effect  of  (i)  the  orbital  vdority  of  the  waves  and  (ii)  wind  speed  have  on  different 
polarized  radar  signals  .  Since  we  did  not  have  a  precise  measurement  of  the  wind  vdodty 
at  the  location  where  data  were  collected,  it  was  not  possible  to  verify  Pidgeon’s  theory. 
More  research  is  needed  to  explain  the  mechanisms  which  produced  this  difference  in 
Doppler  shifts  between  the  two  polarizations. 

We  obtained  similar  results  for  S— band  with  a  reduced  clutter  strength  compared  to 
X-band.  Figure  3a  and  3b  show  the  spectra  for  two  S-band  experiments  for  vertical  and 
horizontal  polarization,  respectively.  These  experiments  were  performed  a  few  minutes 
before  the  X— band  experiments  described  earlier.  The  spectra  were  computed  from  data 
taken  from  the  same  general  areas  as  those  described  in  Figure  2.  Again  we  can  see  that 
the  Doppler  shifts  in  the  sea-clutter  spectrum  for  the  H— POL  data  were  greater  than  those 
for  the  V— POL  data.  In  Section  3.6.3  we  shall  present  some  result  which  could  give  some 
clues  to  this  phenomenon. 

Generally,  we  observed  an  asymmetrical  sea-clutter  spectrum,  with  a  longer 
spectral  tail  in  the  direction  of  the  Doppler  shift.  Symmetrical  spectra  were  observed 
occasionally  in  cross-swell  conditions  where  the  Doppler  shift  was  near  zero.  There  was  no 
noticeable  difference  between  sea-clutter  spectra  for  high-  (15  m)  and  low-resolution  (150 
m)  waveforms. 

Different  spectral  results  were  obtained  for  L— band,  UHF  and  VHF.  Figure  4  shows 
the  spectrum  of  an  L— band,  V— POL  experiment.  There  were  two  distinct  peaks  separated 
by  a  Doppler  frequency  of  about  8  Hz.  The  larger  peak  was  at  about  +6  Hz,  while  the 
other,  about  13  dB  lower,  was  at  -2  Hz. 

In  Figure  5  we  show  the  spectra  of  the  L-band,  V-POL  sea  clutter  at  a  fixed  range 
for  five  look  directions.  The  experimental  parameters  are  shown  in  the  figure.  In  order  to 
show  the  separation  and  the  relative  magnitude  of  the  spectral  peaks  clearly,  the  spectra 
were  offset  by  15  dB  from  one  another.  That  is,  the  magnitude  of  the  spectrum  for  0° 
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SPECTRAL  DENSITY  (dBm2/AHz)  SPECTRAL  DENSITY  (dBm2/AHz) 


DOPPLER  FREQUENCY  (Hz) 


(a)  V-POL 


DOPPLER  FREQUENCY  (Hz) 


(b)  H-POL 

Figure  3.  Typical  low- resolution  S-band  sea-dutter  spectra  in  the  up-swell  and 
near-cross-swell  directions. 
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Figure  4.  A  typical  L-band,  V-POL  low  resolution  sea-clutter  spectrum  observed  in 
the  up-swell  direction. 


look  direction  is  correct  as  appeared  in  Figure  5,  but  the  spectral  magnitudes  for  look 
directions  of  22.5°,  45°,  67.5°  and  90°  should  be  15  dB,  30  dB,  45  dB  and  60  dB  higher, 
respectively,  than  those  shown  in  the  Figure.  The  spectra  have  a  progressively  more 
positive  Doppler  shift  as  the  look  direction  increased  from  0°  to  90°,  indicating  that  the 
up-swell  direction  was  near  90°.  However,  the  separation  between  the  two  peaks  were 
fairly  constant  for  various  look  directions.  Also  in  this  particular  case,  the  spectral  peak 
having  a  greater  positive  Doppler  shift  had  a  greater  magnitude. 

Figure  6  shows  the  spectra  of  an  L-band,  V-POL  sea-clutter  data  set  collected 
from  the  same  general  area  as  those  of  Figure  5  on  a  different  day.  Except  for  the  direction 
of  90°,  the  sea-clutter  spectra  all  had  double  peaks  and  a  negative  mean  Doppler  shift.  In 
this  case,  however,  the  spectral  peak  with  a  greater  negative  Doppler  shift  had  a  greater 
magnitude.  For  the  90°  look  direction,  the  sea-clutter  spectrum  was  nearly  symmetrical, 
indicating  that  this  was  the  cross-swell  direction.  This  condition  was  seldom  observed. 
Usually,  one  of  the  two  peaks  was  substantially  higher  (15  to  20  dB)  than  the  other 
Consequently,  at  longer  ranges  (>  3  km),  the  lower  spectral  peak  would  be  below  noise 
level,  leaving  only  the  higher  spectral  peak  visible. 

For  UHF  and  VHF,  we  also  observed  double  peaks  in  their  spectra,  as  shown  in 
Figure  7  for  a  VHF  experiment.  We  will  briefly  discuss  the  origin  of  double  spectral  peaks 
in  Section  3.4. 
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SPECTRAL  DENSITY  (dBm  /A  HZ) 


DOPPLER  FREQUENCY  (Hz) 

Figure  5  Typical  low- resolution  L-band,  V-POL  sea-clutter  spectra  with  positive 

averaged  Doppler  shifts. 
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SPECTRAL  DENSITY  (dBm  /A  HZ) 


0 


DOPPLER  FREQUENCY  (Hz) 


Figure  6.  Typical  low-resolution  L— band  sea— clutter  spectra  with  negative  averaged 
Doppler  shifts. 
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3.2  Decorrelation  Time  Analysis 

For  a  symmetrical  spectrum,  centred  about  zero  Doppler,  the  autocorrelation 
function  is  a  real  and  even  function.  If  the  spectrum  has  a  net  Doppler  shift,  the 
autocorrelation  function  is  simply  the  product  of  the  real  function  (from  the 
zero-Doppler-shift  case)  with  a  complex  sinusoid,  exp(±jmT),  where  Ud  is  the  angular 
Doppler  frequency  and  T  is  the  time  lag.  The  sign  in  the  complex  sinusoid  is  associated 
with  the  direction  of  the  Doppler  shift.  For  a  positive  Doppler  shift,  the  exponential  is 
exp(-jmT),  and  for  a  negative  Doppler  shift,  the  exponential  is  exp(+ju;dT).  If  one  has  a 
symmetrical  sea-clutter  spectrum  that  can  be  represented  by  a  Gaussian  function,  the 
spectral  width  and  the  Doppler  shift  can  be  determined  from  the  autocorrelation  function. 
Although  in  general,  sea-clutter  spectra  are  seldom  symmetrical,  we  can  obtain  a  fairly 
good  estimate  of  the  sea-clutter  Doppler  shift  and  spectral  width  using  this  method. 

For  X—  and  S— bands,  we  found  the  following  empirical  relationship  works 
reasonably  well:  the  3  dB  spectral  width  of  the  sea-clutter  is  approximately  equal  to  the 
inverse  of  twice  the  decorrelation  time. 

Figures  8a  and  8b  show  the  spectrum  and  the  autocorrelation  function,  respectively, 
of  the  S— band  sea-clutter  from  a  resolution  cell.  In  Figure  8a  we  superimposed  a  smoothed 
curve  (using  a  10— point  block  averaging)  to  enable  us  to  determine  the  Doppler  shift  and 
the  3  dB  spectral  width,  which  we  estimated  to  be  17  Hz  and  23  Hz  respectively.  In  Figure 
8b  we  estimated  the  decorrelation  time  to  be  22  msec.  Using  the  empirical  relationship 
given  in  the  preceding  paragraph,  we  obtain  a  3  dB  spectral  width  of  22.7  Hz.  The  period 
of  the  modulating  sinusoid  for  the  autocorrelation  function  was  determined  from  the 
distance  of  the  second  zero-crossing  of  the  imaginary  part  from  the  origin.  The  imaginary 
part  of  the  autocorrelation  function  started  at  zero,  and  the  point  of  the  second 
zero-crossing  represented  a  full  period.  We  estimated  the  period  from  Figure  8b  to  be 
about  53.5  msec.  This  gave  an  estimated  Doppler  shift  of  18.8  Hz.  These  figures  agreed 
fairly  well  with  the  spectrum  shown  in  Figure  8a. 

Typical  decorrelation  times  for  X-,  S-bands  were  10  msec  and  30  msec, 
respectively,  over  most  of  the  experiments.  Decorrelation  time  analysis  is  less  meaningful 
for  lower  frequencies  owing  to  the  presence  of  double  peaks  in  their  spectra.  In  this  case 
the  autocorrelation  function  provides  information  only  on  the  larger  spectral  peak. 

3.3  Clutter  Coefficient  as  a  Function  of  Grazing  Angle 

Figure  9  shows  the  vs  grazing  angle  profiles  for  the  X— band  data  whose  spectra 
are  compared  in  Figure  2a  and  2b.  The  Doppler  shifts  typical  of  the  sea-clutter  are 
indicated  for  each  look  direction.  For  both  vertical  and  horizontal  polarizations,  a°  was 
substantially  higher  in  directions  of  positive  Doppler  shift.  The  magnitudes  of  a0  were 
similar  for  both  polarizations  in  the  same  look  direction.  The  120°  azimuth  and  the  30° 
azimuth  were  identified  as  the  up-swell  and  cross-swell  directions  based  on  the  estimated 
Doppler  shifts  in  their  respective  spectra.  Generally,  a°  decreased  as  the  look  direction 
moved  away  from  the  up— swell  direction.  We  did  not  observe  further  decrease  in  a0  for 
directions  with  negative  Doppler  shifts  in  the  spectra. 
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CORRELATION  SPECTRAL  DENSITY  (<fBm*/AHz) 


DOPPLER  FREQUENCY  (Hz) 


(a)  Spectrum 


TIME  LAG  (m**c) 

(b)  Autocorrelation  function. 

Figure  8.  Relationship  between  sea-clutter  spectrum  and  autocorrelation. 
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Figure  9.  Typical  <r°  vs.  grazing  angle  profiles  of  X-band  sea  clutter  in  various  look 
directions. 
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We  obtained  similar  results  for  S— band,  as  shown  in  Figure  10.  This  S— band 
experiment  was  performed  approximately  10  minutes  before  those  of  the  X-band 
experiments  described  in  the  preceding  section.  The  S-band  <7°  was  approximately  10  dB 
lower  than  the  corresponding  X— band  values  for  the  same  look  direction.  The 
vertical— to-horizontal  ratio  of  a°  was  also  higher  for  S-band  than  for  the  X— band.  The 
exact  value,  however  was  dependent  on  the  look  direction  with  respect  to  that  of  the  sea 
waves. 


A  different  characteristic  was  observed  for  L— band.  In  Figure  11,  we  show  the  a0 
vs.  grazing  angle  profiles  for  two  L— band  experiments,  one  V— POL  and  one  H— POL.  In 
both  cases,  the  look  directions  were  0°,  22.5°,  45°,  67.5°  and  90°.  The  representative 
Doppler  shift  for  the  sea  clutter  in  each  look  direction  ranged  from  3  Hz  (0°)  to  about  7  Hz 
(90°).  The  o°  vs  grazing  angle  profiles  for  all  five  look  directions  were  very  similar.  For 
horizontal  polarization,  the  values  were  very  low  (about  —70  dB  at  4°  grazing  angle). 
The  anomalous  peaks  in  the  <7°  vs  grazing  angle  profile  were  due  to  targets  of  opportunity 
(boats)  whose  Doppler  fell  within  the  clutter  region,  and  the  rising  trend  for  grazing  angles 
below  1.7°  was  due  to  receiver  noise. 

In  Figure  12  we  show  the  <r°  vs  grazing  angle  profiles  for  an  L— band  experiment 
performed  on  a  different  day.  Here,  the  Doppler  shifts  of  the  spectra  in  all  directions  were 
negative.  Comparing  Figures  11  and  12,  we  found  that  values  of  o°  were  similar  for 
L— band  sea  clutter  with  either  positive  or  negative  Doppler  shifts,  except  for  negative 
Doppler  shifts  near  zero.  This  implied  that  the  up-swell/down-swell  ratio  of  a0  was  close 
to  unity,  and  the  up-swell/cross-swell  ratio  of  a°  was  moderately  high  for  L-band  (about 
10  dB  from  Figure  12).  The  vertical-to-horizontal  ratio  of  o°  was  very  large  at  L-band. 

We  obtained  similar  results  for  UHF  and  VHF.  Figure  13  shows  a  typical  a0  vs 
grazing  angle  profile  for  VHF  sea  clutter.  In  the  UHF  case,  only  vertical-polarization  data 
were  available. 

3.4  Comparison  with  Established  Sea-Clutter  Models 

A  number  of  sea-clutter  models  are  available  to  predict  a0,  and  they  provide  fairly 
reliable  a°  estimates  at  high  grazing  angles  (above  10°).  However,  considerable 
disagreement  exists  among  some  of  the  sea-clutter  models  at  low  grazing  angles  and  at  low 
sea  states.  The  Georgia  Institute  of  Technology  (GIT)  sea-clutter  model[12]  provides 
prediction  of  sea-clutter  from  1  GHz  to  100  GHz.  The  Sittrop  model[13]  is  derived 
from  experimental  data  at  X— band  and  Ku  band.  To  date,  the  data  for  UHF  and  VHF  are 
insufficient  to  develop  comprehensive  models. 

As  discussed  earlier,  the  North  Truro  data  might  include  effects  of  ocean  currents 
which  the  GIT  and  the  Sittrop  models  do  not  address.  Instead  these  models  employ  the 
look  direction  with  respect  to  the  wind  direction  as  an  input  parameter.  There  is  no 
assurance  that  sea  clutter  observed  in  the  presence  of  ocean  current  will  be  similar  to  that 
observed  from  a  wind  driven  sea  surface.  Consequently,  a  direct  comparison  between  the 
North  Truro  results  with  these  models  was  not  possible.  Nevertheless,  it  would  be 
qualitatively  informative  to  compare  the  North  Truro  results  with  predictions  obtained 
from  the  above  mentioned  models  assuming  similar  sea  conditions.  The  reasoning  was 
that,  in  an  open  ocean,  a  sea  wave  travelling  in  a  given  direction  could  be  a  result  of  a 
steady  wind  blowing  in  the  same  direction. 
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Figure  11  Typical  <r°  vs.  grazing  angle  profiles  of  L-band  sea  clutter  in  various  look 
directions  from  up-swell. 
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Figure  12  Typical  <r°  vs.  grazing  angle  profiles  of  L— band  sea  clutter  in  various  look 
directions  from  down-swell. 
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Figure  13  Typical  <r°  vs.  grazing  angle  profiles  for  VHF  sea  clutter  in  various  look 
directions. 


In  Figure  14,  we  plotted  the  North  Truro  X-band  a°  vs  grazing  angle  profiles 
together  with  those  predicted  by  the  GIT  and  the  Sittrop  models.  For  the  GIT  model,  we 
used  an  average  wave  height  of  0.2  m,  and  for  the  Sittrop  model,  we  used  a  sea  state  value 
of  2.  In  the  range  of  grazing  angles  from  1.5°  to  5°,  we  found  a  reasonably  good  agreement 
between  the  North  Truro  X-band  data  (up-swell  condition)  and  the  GIT  model  (upwind 
condition).  The  match  was  less  good  for  grazing  angle  below  1°.  Also  the  GIT  model 
predicted  a  much  smaller  upwind-to-crosswind  c°  ratio  than  the  North  Truro 
up-swell/cross-swell  <r°  ratio.  The  Sittrop  model  yielded  reasonably  close  results  for 
grazing  angles  above  3°  and  below  0.5°.  Also  its  cross-wind  results  matched  the  North 
Truro  cross-swell  results  closer  than  did  the  GIT  model.  However  considerable 
discrepancy  was  found  in  the  range  of  grazing  angle  between  0.5°  and  3°. 

Results  for  S-  and  L— bands  are  shown  in  Figures  15  and  16,  respectively.  Only  the 
GIT  model  was  used  in  comparison. 

Profiles  of  cr°  vs  grazing  angle  for  UHF  and  VHF  sea  clutter  are  given  in  Figure  17. 
These  profiles  were  very  consistent  over  the  two— week  period.  There  were  no  H-POL  data 
for  UHF. 

There  was  a  much  greater  variability  in  <r°  at  low  grazing  angles(<  2°)  than  at 
higher  grazing  angles.  To  show  how  <t°  varied  with  radar  frequency  in  the  low  grazing 
angle  region,  we  plotted  in  Figure  18  the  a0  values  for  various  bands  at  a  grazing  angle  of 
2°.  The  corresponding  V— POL/H— POL  a°  ratio  were  also  shown. 

Examining  Figures  14  to  18,  we  noted  the  following: 

(i)  V-POL  <r°  for  X— ,  S-  and  L-bands  were  similar  except  in  the  up-swell 
direction  where  the  X-band  values  were  about  10  dB  higher; 

(ii)  V— POL/H— POL  ratio  of  <r°  increased  with  decreasing  frequency; 

(iii)  a0  increased  as  frequency  decreased  below  L— band. 

The  above  results  together  with  the  spectral-analysis  results  reinforced  the  theory 
that  the  main  mechanisms  contributing  to  sea  clutter  are  (a)  surface  effects  such  as  water 
sprays  and  droplets,  (b)  wave  facets  and  (c)  Bragg  resonance. 

The  surface  effects  affected  mainly  the  X-band  and  to  a  lesser  extent,  S-band. 
This  could  account  for  the  high  <r°  values  at  X-band  compared  to  those  of  the  S-  and 
L-bands  and  the  relatively  low  V-POL/H-POL  ratio.  Scattering  of  X— band  signals  by 
water  droplet  occurs  in  the  upper  portion  (relative  to  those  of  S—  and  L-bandsl  of  the 
Rayleigh  region  and  will  produce  a  much  stronger  return  than  at  S-  and  L-banas.  The 
spherical  shape  of  water  droplets  means  that  both  V-POL  and  H-POL  signals  were 
scattered  equally  well. 

Rayleigh  scattering  predicts  an  effective  RCS  proportional  to  1/A4  where  A  is  the 
radar  wavelength.  However,  the  measured  <r°  values  for  S-  and  L-bands  did  not  follow 
this  relationship.  Scattering  of  microwave  from  small  wave  facets  of  capillary  waves  may 
account  for  the  less  rapid  rate  of  decrease  of  effective  RCS  with  increasing  wavelength  than 
predicted  by  Rayleigh  scattering  theory. 
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GRAZING  ANGLE  (Deg.) 

Figure  14  Averaged  X-band  vs  grazing  angle  profiles  for  North  Truro  in  the  up-svell  and 

cross-swell  conditions  and  those  calculated  fron  the  GIT  and  Sittrop  nodels  assuring 
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crosswind  conditions. 
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Figure  16  Averaged  L-band  <r°  vs  grazing  angle  profiles  for  North  Truro  in  the  up-swell  and 
cross-swell  conditions  and  those  calculated  froa  the  GIT  aodel  assuming  upwind  and 
crosswind  conditions. 


NORTH  TRURO 
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cross-swell  conditions. 


NORTH  TRURO  SEA  STATE  2  (GRAZING  ANGLE  =  2° ) 
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RADAR  FREQUENCY  (MHz) 

as  a  function  of  frequency  for  the  North  Truro  sea-clutter  data 


A  capillary  wave  is  generated  by  wind,  but  its  velocity  of  propagation  is  controlled 
mainly  by  surface  tension  of  the  water.  The  wavelength  of  a  capillary  wave  is  measured  in 
terms  of  centimeters.  Scattering  from  wave  facets  could  occur  in  the  upper  portion  of  the 
Rayleigh  region  and  the  resonance  region. 

The  double  peaks  exhibited  by  sea  clutter  at  low  frequency  bands  and  the  rising  a° 
below  L— band  suggested  that  a  different  mechanism  was  the  major  contributor  of  sea 
clutter  at  these  frequencies.  This  was  attributed  to  Bragg  resonance  of  surface  waves  [14]. 
The  Bragg  Doppler  shifts  of  sea  surface  waves  at  low  grazing  angles  is  given  approximately 
by  [14]: 


fB  “  ± 


ih 


1/2 


(5) 


where  g  is  the  gravitational  acceleration  (9.8  m/sec2),  and  A  is  the  radar  wavelength. 


Applying  Eqn  (5)  to  the  VHF  case,  we  obtained  a  value  of  ±1.32  Hz,  giving  a  peak 
separation  of  2.64  Hz,  which  was  in  good  agreement  with  the  figure  observed  in  Figure  7. 
Similar  agreements  were  obtained  for  the  L— band  and  UHF  cases. 

3.5  Sea-Clutter  Amplitude  Distribution 

Statistical  properties  of  sea  clutter  are  of  interest  to  radar  designers  for  the 
prediction  of  probability  of  detection(Pd)  and  probability  of  false  alarm(Pfa).  In  this 
sub-section,  results  are  shown  for  the  relative  goodness  of  fit  of  sea-clutter  amplitude 
distributions  to  some  elementary  statistical  models.  We  considered  Rayleigh,  Weibull[15], 
lognormal[16]  and  the  K-distribution[17,18].  The  model  cumulative  distribution  function 
(cdf),  probability  density  function  (pdl),  ana  the  moments  necessary  to  compute  the  model 
parameters  are  defined  in  Equations  (6)  through  (20). 


Ravleieh 

cdf: 

P(x)  =  1  exp{  } 

0  <  x  <  ® 

(6) 

pdf: 

PM  =  T  exp{-  f  } 

(7) 

First 

moment:  <x>  =  ■  — 

(8) 

27 


Weibull 


cdf:  P(x)  =  1  -  exp{-  -f-  } 

(9) 

0  <  X  <  CD 

pdf:  p(x)  =  n[  a  ]  exp{-  fl  } 

(10) 

First  moment:  <x>  =  a1/11  T(1  +  -i~) 

(H) 

second  moment:  <x2>  =  T(1  +  ) 

(12) 

where  T(z)  is  the  gamma  function. 


Lognormal 

cdfc  P(x)  =  1-4-  erfc(  )  (13) 

0  <  X  <  OD 

where  erfc(z)  is  the  complementary  error  function. 


pdf:  p(x)  =  ~h^r, exp{'  JAjLr  m)i  > 

x(za)1M 

(14) 

First  moment:  <x>  =  exp(m  + 

(15) 

Second  moment:  <x2>  =  exp(2m  +  a) 

(16) 

K-Distribution 

cdf:  P(x)  =  1  -  rfjJ  K,(b*) 

0  <  X  <  OD 

(17) 

Pdf-  p(x)  TJVj-  [ — 2 — ]  Ki/-l(bx) 

(18) 
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where  Kj/(z)  is  the  modified  Bessel  function  of  the  third  kind  of  order  v. 


Second 

Fourth 


moment:  <x2> 
moment:  <x*> 


(19) 

(20) 


Note  that  parameters  v  and  b  of  the  K-distribution  are  usually  defined  in  terms  of 
the  first  and  second  moments  of  the  clutter  RCS  which  are  the  second  and  fourth  moments 
of  the  clutter  amplitude,  respectively. 

From  the  radar  perspective,  standard  statistical  tests  such  as  the  chi-square 
goodness  of  fit  test  [19]  are  of  limited  use  for  clutter  data.  The  reason  is  that  these  tests 
place  an  equal  importance  in  the  goodness  of  fit  on  all  regions  of  the  probability  space.  In 
radar  applications,  goodness  of  fit  of  clutter  statistics  to  a  model  is  important  mainly  in  the 
low  Pfa  region. 

In  the  standard  chi-square  goodness  of  fit  test,  the  pdf  of  the  statistical  model  is 
divided  into  K  intervals  of  equal  probability  (i.e.,  1/K).  Parameter  x2  is  defined  as: 


(f»  -  -r>* 


(21) 


where  fi  =  observed  number  of  occurrences  of  a  clutter  sample  having  an  amplitude  within 
the  ith  interval,  and  N  =  total  number  of  amplitude  samples  forming  the  histogram. 

(N/K)  is  the  expected  number  of  occurrences  in  each  interval  for  the  statistical  model. 

A  lower  x2  value  indicates  a  better  fit  to  a  model.  To  determine  the  relative 
goodness  of  fit  of  the  sea  clutter  from  a  group  of  resolution  cells,  we  arbitrarily  assigned  a 
resolution  cell  number  to  each  cell  in  the  group.  The  x2  value  for  data  from  each  cell  was 
computed.  The  results  for  various  models  can  be  plotted  as  a  function  of  the  resolution  cell 
number  to  give  a  visual  indication  of  the  relative  goodness  of  fit  to  the  models. 
Alternatively,  the  ensemble  average  of  x2  for  a  group  of  resolution  cells  can  be  computed 
for  each  model  and  compared  to  give  a  quantitative  measure  of  the  relative  goodness  of  fit. 

We  performed  the  standard  chi-square  goodness  of  fit  test  on  a  set  of  X— band, 
V— POL,  15  m  waveform  resolution  data  using  a  value  of  K=50.  The  length  of  the  time 
series  N  was  10000.  We  obtained  the  following  x2  values:  3190.2,  2637.4,  807.4,  and  7230.6 
for  the  Rayleigh,  Weibull,  lognormal  and  the  K  models,  respectively.  Thus  the  standard 
test  suggested  that  the  data  tor  this  experiment  fitted  the  lognormal  model  best  and  the 
K-distrioution  worst.  Upon  closer  examination,  it  was  found  that  the  data  fitted  the 
K-distribution  poorly  primarily  in  the  Pfa  >  0.1  region.  Figure  19  compares  the  Pfa  vs 
threshold  (Vt)  curve  for  the  above  data  set  with  those  computed  from  the  four  models.  It 
is  seen  that  the  data  fit  the  K-distribution  best  in  the  amplitude  region  for  which  the  Pfa 
is  <  0.1. 
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THRESHOLD  (LINEAR) 


Figure  19  Comparison  of  the  Pfa  characteristics  of  an  X-band,  V-POL  sea-clutter 
data  set  with  those  predicted  by  the  Rayleigh,  Weibull,  lognormal  and  K 
models. 


Modified  Chi-Sauare  Test 

To  obtain  a  procedure  which  brings  out  the  relative  merits  of  various  models  in  the 
important  low  Pfa  region,  we  modified  the  chi-square  test  as  follows.  Referring  to  Figure 
20,  boundaries  of  equal— probability  intervals  were  determined  for  each  of  the  statistical 
models  in  the  amplitude  region  for  which  Pfa  is  <  0.1.  These  boundaries  divided  this  region 
into  K  equal-probability  intervals,  each  with  a  probability  of  0.1/K.  We  defined  the 
modified  x2  index: 
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(22) 


,  v  [fi  •  C-V-)]1 

*• =  k  ~V^~) 


Figure  20  Illustration  of  the  division  of  equal— probability  intervals  for  the  modified 
chi-square  test. 


This  is  equivalent  to  assigning  a  zero  weighting  in  the  amplitude  region  where  the  Pfa  is 
greater  than  0.1  and  a  uniform  weighting  in  each  of  the  K  intervals.  The  modified 
chi-square  test  is  simply  a  comparison  of  the  values  calculated  for  each  of  the  statistical 

models. 

Using  the  modified  chi-square  test,  goodness  of  fit  comparisons  were  made  among 
sea-clutter  data  for  various  combinations  of  frequency  bands  (X-,  S-  and  L-band), 
polarization  (V-  and  H— POL)  and  waveform  resolutions  (15m  ana  150  m).  For  L— band, 
H-POL  data,  the  sea-clutter  component  was  very  weak,  and  the  composite  returns  were 
dominated  by  receiver  noise.  Consequently,  L-band,  H-POL  sea-clutter  presented  a  good 
fit  to  the  Rayleigh  and  Weibull  models. 
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Figure  21  shows  a  typical  result  of  the  modified  chi-square  test  for  low  resolution 
vertical  polarization  data  (X— band,  V-POL,  150m  waveform).  In  this  figure,  the  modified 
chi-square  indices  for  40  resolution  cells  are  shown.  Cell  "0"  was  at  a  range  of  1  km,  and 
each  point  in  the  figure  represents  a  range  increment  of  150  m.  There  was  a  number  of 
resolution  cells  with  clutter  data  having  a  poor  fit  to  the  Rayleigh  model,  suggesting  that 
non— Rayleigh  sea-clutter  statistics  were  common  even  at  low  waveform  resolutions.  On 
average,  these  data  had  the  best  fit  to  the  Weibull  and  K-distributions  and  the  worst  fit  to 
the  lognormal  model. 

In  Figure  22,  We  plotted  the  Pfa  vs  Vt  curves  for  data  in  cell  No.  23  of  Figure  21. 
The  modified  chi-square  test  showed  that  this  data  set  had  a  poor  fit  to  the  Rayleigh 
model  and  a  reasonably  good  fit  to  the  lognormal  model.  This  was  confirmed  by  the  results 
in  Figure  22,  and  it  also  showed  that  the  data  fitted  the  K-distribution  best. 

A  typical  result  of  the  modified  chi-square  test  for  high  resolution  horizontal 
polarization  data  is  given  in  Figure  23  (S-band,  H— POL,  15  m  waveform).  Cell  "0"  in  this 
figure  was  2  km,  and  each  point  represents  a  range  increment  of  15  m.  In  this  case,  the 
modified  chi-square  indices  for  all  resolution  cells  assuming  a  Rayleigh  model  were  very 
large  (implying  a  poor  Rayleigh  fit)  and  were  beyond  the  scale  employed  in  Figure  23.  The 
data,  however,  had  a  reasonably  good  fit  to  the  lognormal  distribution  in  all  resolution  cells 
included  in  the  test. 

Figure  24  compares  the  Pfa  vs  Vt  curve  for  data  from  a  typical  resolution  cell  with 
those  computed  from  the  models.  The  data  fitted  the  K-distribution  best. 

To  establish  the  relative  merits  of  the  statistical  models  considered  in  this  study,  we 
computed  the  ensemble  average  of  the  modified  chi-square  index  over  a  group  of  resolution 
cells.  We  tested  data  from  various  combinations  of  frequency,  polarization  and  waveform 
resolution.  The  results  are  tabulated  in  Table  II. 

In  nearly  all  cases,  we  obtained  the  best  fit  of  the  data  to  the  K— distribution. 
Generally,  the  fit  to  a  lognormal  distribution  improved  for  high  resolution  (15m) 
waveform  and  for  horizontal  polarization.  The  fit  to  a  Rayleigh  and  Weibull  distributions 
improved  for  low— resolution  waveform  (150  m)  and  vertical  polarization. 

The  results  presented  in  Table  II  together  with  those  of  Figures  21  to  24  indicated 
the  following: 

(i)  The  K-distribution  is  the  most  appropriate  model  for  sea  clutter  in  the 
low  Pfa  region; 

(ii)  Sea-clutter  statistics  approached  those  of  a  Rayleigh  model  for  low 
resolution  and  vertically  polarized  waveforms; 

(iii)  Non-Rayleigh  sea-clutter  statistics  were  common  even  for  relatively  low 
resolution  and  vertical  polarization  waveforms; 

(iv)  The  Weibull  distribution  was  also  a  good  model  for  low-resolution  sea 
clutter. 

(v)  Sea-clutter  statistics  approached  those  of  a  lognormal  model  for  high 
resolution  and  horizontal  polarization  waveforms. 
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PROBABILITY  OF  FALSE  ALARM 


THRESHOLD  (LINEAR) 


Figure  22  Comparison  of  the  Pfa  characteristics  of  a  typical  low  resolution  X-band, 
V— POL  sea-clutter  data  set  with  those  predicted  by  the  Rayleigh,  Weibull, 
lognormal  and  K  models. 
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Figure  23  Comparison  of  the  modified  chi-square  index  of  S-band,  H-POL,  15  m  waveform  sea  clutte 
computed  for  the  Rayleigh,  Veibull,  lognormal  and  K  models  for  a  group  of  resolution 
cells. 


THRESHOLD  (LINEAR) 


Figure  24  Comparison  of  the  Pfa  characteristics  of  a  typical  high  resolution  S-band, 
H-POL  sea-clutter  data  set  with  those  predicted  by  the  Rayleigh,  Weibull, 
lognormal  and  K  models. 
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Table  II.  Summary  of  modified  chi-square  test  on  sea  clutter 


Averaged  nodified  index 

Band 

POL 

Bes. 

(m) 

Rayleigh 

Veibull 

Lognormal 

K 

Sample 

size 

V 

150 

18.11 

10.57 

22.96 

9.46 

80 

X 

H 

150 

168.77 

15.45 

11.68 

8.93 

80 

ms 

15 

150.2 

9.99 

15-18 

7.51 

160 

H 

15 

802.7 

24.58 

12.50 

14.62 

160 

n 

150 

23.01 

18.30 

36.33 

17.24 

80 

s 

H 

150 

361.80 

20.52 

19.14 

13.00 

80 

n 

15 

85.60 

8.09 

18.85 

7.85 

160 

H 

15 

423.01 

31.19 

23.38 

18.99 

160 

T 

V 

150 

15.43 

12.02 

23.42 

11.85 

80 

Li 

V 

15 

19.57 

16.89 

30.49 

17.03 

80 

3.6  Sea  Swell  and  Sea  Spikes 

One  of  the  strange  phenomena  observed  in  sea  clutter  is  that  if  sea  spikes  —  random 
occurrences  of  abrupt  increase  of  sea-clutter  magnitude  of  time  duration  in  the  order  of  a 
second.  Ward  [18]  identifies  a  component  in  the  sea  clutter  having  a  spatially  varying 
mean  level  as  the  result  of  the  bunching  of  scatterers  by  sea-swell  structures.  Thus  sea 
swell  could  also  be  responsible  for  the  occurrence  of  sea  spikes.  The  verification  of  this 
theory  requires  that  the  sea-clutter  data  be  analyzed  together  with  ground  truth  data 
(regarding  the  wave  crest  and  heights,  etc.)  taken  simultaneously.  This  type  of  ground 
truth  data  was  not  available  for  the  North  Truro  experiment.  Nevertheless,  there  were 
some  high  resolution  data  suitable  for  examining  sea-spike  phenomena. 

3.6.1  Range-time-intensity  plots. 

A  useful  tool  in  analyzing  sea  clutter  both  spatially  and  temporally  (in  a  time  scale 
measured  in  seconds')  is  the  range-time-intensity  (RTI)  plot.  Here  high-resolution  f  15  m) 
sea-clutter  data  taken  from  a  number  of  resolution  cells  in  a  given  direction  tor  an 
extended  period  of  time  (20  seconds)  were  used  to  form  a  RTI  plot. 
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The  sample  mean  of  the  data  was  removed  to  suppress  any  possible  ground  clutter 
contribution  from  the  antenna  sidelobes.  The  squared  magnitudes  of  the  resulting  time 
series  were  then  summed  over  contiguous  blocks  of  100  samples.  The  result  is  a  100— point 
time  series  with  each  point  representing  the  averaged  clutter  magnitude  over  a  200  msec 
time  interval.  This  is  done  for  data  of  100  contiguous  resolution  cells.  The  resulting  data 
set  is  then  plotted  in  a  3-D  format. 

3.6.2  Vertically  and  horizontally  polarized  sea  clutter  at  high  resolution. 

A  typical  RTI  plot  for  an  X— band,  V-POL  sea-clutter  data  set  is  shown  in  Figure 
25.  In  this  figure  the  x— dimension  is  range  (100  points  —  15  m  per  point),  the  y-dimension 
is  time  (100  points  —  200  msec  per  point)  and  the  z-dimension  is  the  equivalent  radar  cross 
section  (mJ). 

Spectral  analysis  of  this  data  set  showed  that  the  sea-clutter  spectrum  had  a  mean 
Doppler  shift  of  50  Hz,  indicating  that  the  water  waves  were  moving  towards  the  radar  at  a 
speed  of  about  0.82  m/sec.  The  RTI  plot  of  this  data  set  clearly  showed  the  regular 
patterns  associated  with  the  wave  crests  and  troughs.  If  we  start  from  a  certain  resolution 
cell  with  a  relatively  high  clutter  magnitude  a  time  zero,  we  can  detect  a  gradual  migration 
of  the  high  clutter  magnitude  into  progressively  closer  resolution  cells. 

From  Figure  25,  we  measured  approximately  the  distance  over  which  a  given  high 
clutter  magnitude  migrated  in  20  seconds  to  be  234  m.  This  corresponded  to  a  phase 
velocity  of  11.7  m/sec  for  the  water  wave.  This  phase  velocity  was  not  detected  from 
spectral  analysis  because  water  particles  or  wave  facets  do  not  travel  at  the  phase  velocity. 
Indeed  a  velocity  of  11.7  m/sec  should  give  a  Doppler  shift  of  about  780  Hz  at  X-band. 
After  accounting  for  aliasing,  it  did  not  match  that  obtained  from  spectral  analysis  (about 
+50  Hz).  The  Doppler  shift  measured  from  the  sea-clutter  spectrum  represents  the 
velocity  of  the  water  particles  or  wave  facets,  probably  resulting  from  the  effect  of  ocean 
currents. 

Figure  26  shows  the  RTI  plot  for  an  X-band,  H— POL  data  set.  This  experiment 
was  performed  about  6  minutes  after  the  X-band,  V-POL  experiment  (Figure  25).  In  this 
case  we  did  not  observed  the  regular  wave-like  pattern  as  we  did  in  the  case  of  the  vertical 
polarization  case. 

The  large  clutter  magnitudes  of  relatively  short  duration  observed  in  isolated 
resolution  cells  were  associated  with  sea  spikes.  In  general,  the  H— POL  sea  spikes  were 
significantly  higher  than  the  V-POL  sea  spikes.  However,  the  magnitudes  of  the  sea 
clutter  during  the  time  intervals  between  two  sea  spikes  were  much  lower  for  horizontal 
polarization  compared  with  those  of  the  vertical  polarization.  In  Section  3.6.3  we  shall 
show  the  time  waveforms  of  sea-clutter  time  series  containing  sea  spikes  and  examine  in 
greater  details  the  spectrum  of  the  sea  spikes. 

Figure  27  shows  the  RTI  plot  of  another  X-band,  V-POL  data  set.  In  this  case  the 
look  direction  was  perpendicular  to  that  for  the  data  set  depicted  in  Figure  25.  Since  the 
look  direction  was  near  cross-swell,  the  regular  wave-like  patterns  were  not  visible. 
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Figure  25  A  typical  RTI  plot  of  X-band,  V-POL,  15  m  resolution  sea-clutter  data  observed  in  the 
up-svell  direction. 
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Figure  26  A  typical  RTI  plot  of  X-band,  H-P01,  15  m  resolution  sea-clutter  data  observed  in  the 
up-svell  direction. 


X-BAND,  V-POL,  15  m  AZ  ■  30  0.025 


d  (zw)  A1ISN31NI  o 


41 


Figure  27  A  typical  RTI  plot  of  X-band,  V-POL,  15  n  resdiution  sea-clutter  data  observed  in  the 
cross-swell  direction. 


Similar  results  were  obtained  for  the  S— band  cases.  Figures  28  and  29  show  the  RTI 
plots  for  an  S-band,  V-POL  and  a  H-POL  sea-clutter  data  set,  respectively.  These 
experiments  were  performed  about  10  minutes  before  the  X-band  experiments  (Figures  25  - 
27).  Again  the  regular  wave-like  pattern  was  observed  in  the  V— POL  case,  while  a  more 
spiky  and  random  pattern  was  observed  in  the  H— POL  case.  The  phase  velocity  of  the  sea 
swell  estimated  from  the  S— band,  V-POL  RTI  plot  was  in  good  agreement  with  those 
estimated  from  the  corresponding  X— band,  V-POL  RTI  plots. 

3.6.3  Temporal  and  Spectral  Characteristics  of  Sea  Spikes. 

In  this  section  we  show  several  sea-clutter  time  series  containing  sea  spikes.  Figure 
30  shows  the  in-phase  waveform  of  an  X-band,  V-POL,  15  m  resolution  sea-clutter  time 
series,  together  with  the  corresponding  spectra  computed  for  each  two-second  interval 
(called  sectional  spectra).  The  experiment  parameters  are  shown  in  the  figure.  We 
observed  several  time  intervals  of  duration  in  the  order  of  two  seconds  in  which  the 
sea-clutter  amplitude  increased  appreciably  (e.g.,  pulses  No.  3000  to  3800;  5000  to  6200; 
7300  to  8000).  The  Doppler  shifts  and  the  spectral  spread  of  the  spectra  computed  for 
these  segments  of  waveforms  were  substantially  greater  compared  with  those  computed  for 
segments  with  a  lower  clutter  amplitude. 

Figure  31  shows  the  in— phase  waveform  and  the  corresponding  sectional  spectra  for 
an  X— band,  H— POL,  15  m  resolution  sea-clutter  time  series.  In  this  case  the  sea-clutter 
amplitudes  were  fairly  low  during  the  20-second  time  interval,  except  between  pulses 
No.3100  to  3400,  where  a  large  spike  appeared.  The  maximum  instantaneous  amplitude  of 
the  sea  spike  was  almost  3  m,  which  corresponded  to  an  RCS  of  about  9  m2.  We  observed 
that  the  power  and  the  Doppler  shift  of  the  spectrum  for  the  waveform  segment  containing 
the  sea  spike  was  significantly  greater  than  those  of  the  other  sectional  spectra. 

The  results  of  Figures  30  and  31  suggested  the  following: 

(a) .  The  scatterers  responsible  for  the  occurrence  of  sea  spikes  have  a  greater 

velocity  than  those  for  ordinary  sea  clutter. 

(b) .  The  spectral  power  of  H-POL  sea  clutter  was  dominated  by  that  of  the  sea 

spikes. 

Item  (a)  lends  support  to  the  theory  associating  the  occurrence  of  sea  spikes  with 
sea  swells.  The  sea  waves  will  likely  experience  break-up  at  the  crest  of  a  sea  swell.  The 
large  number  of  water  particles  thrown  up  by  the  breaking  waves  could  explain  both  the 
greater  magnitude  and  Doppler  shifts  (at  X—  and  S— bands)  of  sea  spikes. 

Item  (b)  could  explain  the  observed  greater  Doppler  shift  of  the  horizontal  sea 
clutter  compared  with  that  in  vertical  polarization.  If  the  contribution  of  H— POL  spectral 
power  comes  mainly  from  those  of  sea  spikes,  then  the  average  Doppler  shifts  of  a  H— POL 
sea-clutter  spectrum  will  be  greater. 

Figures  32a  and  32b  show  the  averaged  spectrum  of  the  sea-clutter  time  series 
depicted  in  Figures  30  and  31,  respectively.  It  is  clear  that,  while  the  Doppler  shifts  of  the 
sea  spikes  were  similar  for  both  vertical  and  horizontal  polarization,  the  averaged  Doppler 
shift  of  the  spectrum  in  Figure  32b  (H-POL)  is  greater  than  that  of  the  spectrum  in  Figure 
31a  (V-POL). 
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Figure  28  A  typical  RTI  plot  of  S-band,  V-POL,  15  ■  resolution  sea-clutter  data  observed  in  the 
up-svell  direction. 
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Figure  29  A  typical  RTI  plot  of  S-band,  H-POL,  15  m  resoution  sea-clutter  data  observed  in  the 
up-svell  direction. 
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data  set  containing  sea  spikes. 
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data  set  containing  sea  spikes. 


4. 


SUMMARY  AND  CONCLUSIONS 


An  extensive  analysis  was  performed  on  sea-clutter  data  collected  in  North  Truro, 
Massachusetts.  These  data  facilitated  the  examination  of  sea  clutter  over  a  wide  range  of 
frequencies  and  waveform  resolutions  under  similar  sea  conditions.  The  results  of  the 
analysis  supported  the  established  theory  that  the  three  main  mechanisms  contributing  to 
sea  clutter  are  water  droplets(sea  spray),  wave  facets  and  Bragg  resonance.  The  analysis 
provided  insights  into  the  conditions  under  which  these  mechanisms  become  dominant. 
For  example,  contribution  to  sea  clutter  of  scattering  by  water  sprays  and  droplets  is 
important  at  high  frequency  bands  (X  and  S).  At  L— band  and  lower  the  Bragg  resonance 
effect  began  to  dominate  as  evidenced  by  the  concentration  of  spectral  power  around  the 
characteristic  Bragg  Doppler  shifts. 

We  compared  the  <r°  vs  grazing  angle  results  to  the  GIT  and  Sittrop  sea-clutter 
models.  Due  to  the  limitations  in  the  ground-truth  data,  we  did  not  draw  any  conclusions 
with  regard  to  the  relative  merits  of  these  models.  Instead,  we  presented  the  results 
together  with  some  results  calculated  from  the  two  models  so  that  they  can  be  qualitatively 
compared.  More  work  is  required  to  validate  these  predictive  sea-dutter  models  in  the 
low— angle  region,  particularly  in  conditions  where  ocean  current  is  present. 

Results  of  statistical  analysis  showed  that  sea-clutter  amplitudes  from  a  fixed 
resolution  cell  fit  a  K-distribution  for  all  combinations  of  frequencies  (among  X-,  S-  and 
L— bands),  polarization  (V  and  H)  and  waveform  resolutions  (15m  and  150m).  The 
K-distrif>ution  fitting  the  data  approached  the  limiting  cases  of  a  Rayleigh  or  a  lognormal 
model  depending  on  polarization,  waveform  resolution  and  frequency.  For  vertically 
polarized,  low— resolution  waveforms,  sea-clutter  statistics  approached  those  of  a  Rayleigh 
model.  For  horizontal  polarization  and  high  resolution  waveforms,  sea-clutter  statistics 
approached  those  of  a  lognormal  distribution. 

Spatial  and  temporal  characteristics  of  sea  spikes  were  examined  through  the  use  of 
RTI  plots  of  high  resolution  sea  clutter.  The  results  supported  the  theory  that  sea  spikes 
are  associated  with  the  wave  crests  in  sea  swells.  This  association  was  more  evident  from 
the  V— POL  RTI  plots.  For  high  resolution  H-POL  sea  clutter,  the  clutter  magnitudes 
were  substantially  lower  than  those  observed  with  a  V— POL  radar;  however,  the 
magnitude  of  the  sea  spikes  was  generally  higher  than  those  in  the  V— POL  case. 

The  theory  explaining  the  differences  between  the  V-POL  and  H-POL  sea  spikes  is 
still  incomplete.  More  data  with  better  ground  truth  are  needed  to  develop  the  theory  and 
to  improve  sea-clutter  models.  One  possible  way  to  exploit  this  knowledge  of  sea-spike 
characteristics  is  to  use  a  high-resolution  radar  with  polarization  diversity.  Radar  pulses 
are  transmitted  alternately  with  vertical  and  horizontal  polarization.  The 
vertical— polarization  sea  clutter  can  be  used  to  obtain  the  sea-swell  pattern  which  in  turn 
can  be  used  to  predict  the  time  of  arrival  of  a  wave  crest.  Target  detection  will  be 
performed  using  the  horizontally  polarized  returns.  The  detection  threshold  settings  can  be 
adaptively  varied  according  to  the  estimated  time  of  arrival  of  the  wave  crests. 

The  results  presented  in  this  report  should  be  of  use  for  maritime  surveillance  and 
shipboard  radars  where  the  grazing  angle  is  generally  low.  However,  it  should  be 
emphasized  that  these  results  were  observed  at  relatively  low  sea  states  and  should  not  be 
extrapolated  to  higher  sea  states  without  verification.  Data  at  low  grazing  angles  are 
scarce,  particularly  at  lower  frequencies. 
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The  North  Truro  sea-clutter  data  which  were  coherent  and  collected  in  long-dwell 
mode  provided  reliable  a°  estimates  for  cases  of  very  low  clutter  magnitudes  (L— band, 
H— POL).  This  should  provide  more  accurate  information  for  low-angle  sea-clutter 
modeling  than  available  at  present. 

In  conclusion,  through  the  analysis  of  the  North  Truro  data,  we  have  obtained 
valuable  information  and  additional  insights  into  sea-clutter  behaviour  at  low  sea  state 
and  low  grazing  angle.  However  more  data  are  needed,  particularly  at  higher  sea  states,  to 
develop  and  verify  sea-clutter  models. 
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1 3.  ABSTRACT  |  ■  brief  end  tactual  summary  of  the  document  It  may  also  wear  elsewhere  in  the  body  o<  the  document  itself,  it  is  highly 
desirable  that  the  abstract  of  classified  documents  be  unclassified.  Each  paragraph  of  the  abstract  shall  begin  with  an  indication  of  the 
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In  this  report  ve  present  the  results  of  the  analysis  of  sea-clutter  data  collected  at  North 
Truro,  Massachusetts  by  the  IIT  Lincoln  Laboratory.  The  data  were  obtained  in  conditions  typical  of 
sea  state  2  at  I-,  S-,  L-bands,  CHF  and  VHF.  Ve  exaained  the  data  both  spectrally  and  statistically. 
The  results  shoved  that  sea  clutter  exhibited  very  different  spectral  characteristics  at  higher 
frequencies  compared  to  those  at  lover  frequencies.  Experimental  sea-clutter  coefficients  as  a 
function  of  grazing  angle  at  various  frequency  bands  vere  obtained  for  up-svell  and  cross-swell 
conditions,  lost  sea-clutter  models  use  the  radar  look  direction  relative  to  the  vind  direction  as  an 
input  parameter  fe.g. ,  up-vind,  cross-wind,  etc) .  Hence  a  direct  comparison  of  sea-clutter  magnitude 
vith  those  calculated  from  aodels  under  different  vind  conditions  vas  not  made.  Instead,  ve  showed 
these  results  together  vith  those  calculated  from  two  -  sea-clutter  models,  assuming  upvind  and 
cross-vind  conditions.  The  data  vere  quite  dissimilar  in  the  lov-grazing-angle  regions  to  those 
predicted  by  the  tvo  aodels.  Further  research  on  sea  clutter  in  the  lov-grazing-angle  regions  is 
needed  to  verify  and  explain  the  observed  differences.  Ve  found  the  sea-clutter  amplitude  statistics 
fit  the  I-distribution  best  in  the  important  lov-probability-of-false-alarm  region.  The 
distribution  approached  that  of  a  Rayleigh  model  for  vertical  polarization  and  lov-resolution 
waveforms.  However,  there  was  a  significant  proportion  of  resolution  ceils  with  sea  clutter  having  non-Rayleigh 
statistics  even  in  relatively  low— resolution  waveforms  and  vertical  polarization.  For  horizontal 
polarization  and/or  high  resolution  waveforms,  the  statistics  approached  those  of  a  lognormal  model. 
The  statistical  results  provided  further  evidence  that  the  K-distribution  can  serve  as  a  limiting 
distribution  for  sea  clutter.  Spatial  and  temporal  characteristics  of  sea  spikes  vere  examined 
through  the  use  of  range— time— intensity  plots  of  high— resolution  sea  clutter.  The  results  supported 
the  theory  that  sea  spikes  are  associated  vith  the  wave  crests  in  sea  swells.  This  association  was  more 
evident  from  high-resolution,  V-POL  sea-clutter  data.  The  results  of  this  report  are  pertinent  to 
shipboard  radar  applications  where  the  grazing  angle  is  low.  They  also  provide  additional  information 
for  sea-clutter  modeling  in  the  low-angle  region,  particularly  at  the  lower  frequency  bands. 
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